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Vortical structures of turbulence, i.e., vortex tubes and sheets, are studied using one-dimensional 
velocity data obtained in laboratory experiments for duct flows and boundary layers at microscale 
Reynolds numbers from 332 to 1934. We study the mean velocity profile of intense vortical struc- 
tures. The contribution from vortex tubes is dominant. The radius scales with the Kolmogorov 
length. The circulation velocity scales with the rms velocity fluctuation. We also study the spatial 
distribution of intense vortical structures. The distribution is self-similar over small scales and is 
random over large scales. Since these features are independent of the microscale Reynolds number 
and of the configuration for turbulence production, they appear to be universal. 
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I. INTRODUCTION 

Turbulence contains various classes of structures that 
are embedded in the background random fluctuation. 
They are important to intermittency as well as mixing 
and diffusion. Of particular interest are small-scale struc- 
tures, which could have universal features that are inde- 
pendent of the Reynolds number and of the large-scale 
flow. We explore such universality using velocity data 
obtained in laboratory experiments. 

We focus on vortical structures, i.e., vortex tubes and 
sheets. The former is often regarded as the elementary 
structure of turbulence p^ 2 -^ At low microscale Reynolds 
numbers, Re> < 200, direct numerical simulations de- 
rived basic parameters of vortex tubes £i2&£aUL The radii 
are of the order of the Kolmogorov length 77. The total 
lengths are of the order of the correlation length L. The 
circulation velocities are of the order of the rms velocity 
fluctuation (u 2 ) 1 / 2 or the Kolmogorov velocity uk- Here 
(•) denotes an average. The lifetimes are of the order of 
the turnover time for energy-containing eddies L/ [u 2 ) 1 / 2 . 

For these vortical structures, however, universality 
has not been established because the behavior at high 
Reynolds numbers has not been known. At Re> ^ 200, 
a direct numerical simulation is not easy for now. 

The promising approach is velocimetry in laboratory 
experiments. A probe suspended in the flow is used 
to obtain a one-dimensional cut of the velocity field. 
The velocity variation is intense at the positions of in- 
tense structures. Especially at the positions of intense 
vortical structures, the variation of the velocity compo- 
nent that is perpendicular to the one-dimensional cut is 
intense. 9,10 Thus, the velocity variation offers some infor- 
mation about intense structures, although it is difficult 
to specify their geometry. 

The above approach was taken in several 
studiesj 11 ! 12 ' 13 ! 14 ' 15 For example, using grid turbulence 14 
at ReA = 105-329 and boundary layers 1 ^ at Re> = 295- 
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1258, we studied the scale, intensity, and spatial 
distribution of vortical structures. However, the 
Reynolds number could be increased still more. The 
dependence of those features on the large-scale flow, i.e., 
on the configuration for turbulence production, has not 
been known. 

We accordingly use duct flows and boundary layers 
to compare features of intense vortical structures along 
the Reynolds number. The highest Reynolds number 
is ReA = 1934, which exceeds those of the prior stud- 
ies. Our experiments are described in Sec. [Til We dis- 
cuss in Sec. IIIII what information is available from one- 
dimensional velocity data. The scale, intensity, and spa- 
tial distribution of vortical structures are studied in Sees. 
IIVIIVI1 The conclusions are summarized in Sec. IVIII 



II. EXPERIMENT 
A. Duct flow 

The experiment was done in a wind tunnel of the Mete- 
orological Research Institute. Table UJ lists experimental 
conditions and turbulence parameters. The microscale 
Reynolds number ranges from Rca = 719 to 1934. 

We use coordinates x, y, and z in the streamwise, 
spanwise, and floor- normal directions. The correspond- 
ing flow velocities are U + u, v, and w. Here U is the 
mean while u, v, and w are the fluctuations. The origin 
x = y = z = m is on the tunnel floor at the entrance 
to the test section. Its size was Sx — 18 m, Sy = 3 m, 
and Sz = 2 m. We placed a rectangular duct with width 
Sy = 1.3 m and Sz = 1.4 m at x = —2 m. The duct center 
was on the tunnel axis. 

For the flow velocities U* from 11 to 55ms _1 at the 
duct exit, we measured the streamwise and spanwise ve- 
locities at x = 15.5 m and z = 0.6 m. The flow was 
turbulent there. During each of the measurements, the 
flow temperature was constant within ±1°C. 

We used a hot-wire anemometer, which was composed 
of a constant-temperature system and a crossed-wire 
probe. The wires were made of tungsten, 5 /im in di- 
ameter, 1.25mm in sensing length, 1.4mm in separation, 
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TABLE I: Experimental conditions and turbulence parameters: duct-exit or incoming-flow velocity U*, coordinates x and 
z of the measurement position, mean streamwise velocity U, sampling frequency f s , kinematic viscosity u, mean energy 
dissipation rate (e) = 15u((d x v) 2 ) /2, rms velocity fluctuations (u 2 ) 1 ^ 2 and (v 2 ) 1 ^ 2 , Kolmogorov velocity uk = (^(e)) 1 ^ 4 , rms 
spanwise-velocity increment over the sampling interval {Sv 2 ) 1 ^ 2 = {[v(x + U/2f a ) — v(x — U/2f s )] 2 ) 1 ^ 2 , correlation lengths 
L u = f£°(u(x + r)u(x))/(u 2 )dr and L v = / Q °° (v(x + r)v(x)) / (v 2 )dr, Taylor microscale A = [2{v 2 )/ {(d x v) 2 )] 1/2 , Kolmogorov 
length rj = (^ 3 /{e)) 1//4 , and microscale Reynolds number Re^ = Xlv 2 ) 1 ^ 2 /v. The velocity derivative was obtained as d x v = 
[8v(x + r) - 8v(x - r) - v(x + 2r) + v(x - 2r)]/12r with r = U/f s . 
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12.5 


12.5 
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12.5 


12.5 


12.5 
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m 


0.60 


0.60 


0.60 


0.60 


0.60 


0.35 


0.35 


0.30 


0.30 


0.25 


0.25 


U 


m s -1 


4.26 


8.66 


13.0 


17.3 


21.2 


1.51 


3.12 


5.84 


8.75 


10.8 


13.6 


fs 


kHz 


12 


34 


60 


80 


100 


4 


10 


22 


38 


60 


74 


V 


cm 2 s _1 


0.141 


0.141 


0.142 


0.142 


0.142 


0.141 


0.142 


0.142 


0.142 


0.142 


0.143 


(e) 


m 2 s- 3 


0.405 


2.87 


9.13 


19.2 


34.4 


0.0332 


0.226 


1.85 


5.37 


13.8 


24.6 


(u 2 ) 1 ' 2 


m s~ 


0.694 


1.38 


2.11 


2.84 


3.46 


0.283 


0.582 


1.18 


1.80 


2.46 


3.14 


(v 2 ) 1 ' 2 


m s" 1 


0.666 


1.34 


2.04 


2.69 


3.32 


0.242 


0.475 


0.973 


1.46 


1.98 


2.51 


UK 


m s _1 


0.0489 


0.0798 


0.107 


0.128 


0.149 


0.0262 


0.0423 


0.0716 


0.0934 


0.118 


0.137 


(Sv 2 ) 1 / 2 


m s~ 


0.0219 


0.0421 


0.0637 


0.0918 


0.121 


0.00667 


0.0143 


0.0349 


0.0517 


0.0646 


0.0882 
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cm 


55.4 


55.8 


51.9 


46.3 


42.2 


43.9 


48.5 


41.5 


43.4 


47.3 
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L v 


cm 


14.9 


14.5 


14.7 


14.2 


14.2 


6.67 


6.87 


6.33 


5.83 


5.71 


6.11 
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cm 


1.52 


1.15 


0.986 


0.895 


0.826 


1.93 


1.46 


1.04 


0.919 


0.776 


0.742 
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cm 


0.0288 


0.0177 


0.0133 


0.0111 


0.00955 


0.0539 


0.0335 


0.0198 


0.0152 


0.0120 
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Re A 




719 


1098 


1416 


1693 


1934 


332 


488 


716 


945 


1080 


1304 



and oriented at ±45° to the streamwise direction. 

The signal was linearized, low-pass filtered at 24 dB per 
octave, and then digitally sampled at 16-bit resolution. 
We set the sampling frequency / s as high as possible, on 
condition that high-frequency noise was not significant 
in the power spectrum. The filter cutoff frequency was 
one-half of the sampling frequency. For f s < 50 kHz, 
we obtained 10 s data of the streamwise and spanwise 
velocities. For / s > 50 kHz, we obtained 4 x 10 8 data 
of the spanwise velocity alone. This is due to a f s limit 
of our sampling device. Supplementarily, we obtained 
4 x 10 7 data of the streamwise and spanwise velocities at 
f s = 50 kHz. 

The temporal variations were converted into the spa- 
tial variations using Taylor's frozen-eddy hypothesis, 
which requires that the turbulence strength (it 2 ) 1 / 2 /U is 
small enough. This requirement is satisfied in our experi- 
ments where (u 2 ) 1 / 2 /U < 0.2. Even for (u 2 ) 1 / 2 /U ~ 0.3, 
the validity of Taylor's hypothesis was confirmed using 
data obtained simultaneously with two probes separated 
by streamwise distances^ 

By assuming local isotropy ((d x v) 2 ) ~ 2((d x u) 2 ), the 
small-scale statistics were obtained from the spanwise- 
velocity derivative d x v, instead of the usual streamwisc- 
velocity derivative d x u. At small scales, the u component 
measured by a crossed-wire probe is contaminated with 
the w component that is perpendicular to the two wires 
of the probefi 7 - The v component is free from such con- 
tamination. 



B. Boundary layer 

The experiment was done in the same wind tunnel 
with the same instruments as for the duct flows. Table U 
lists experimental conditions and turbulence parameters. 
They are close to those in our prior study^ except that 
the sampling frequency is higher and the data are longer. 
The microscale Reynolds number ranges from Re> = 332 
to 1304. 

Over the entire floor of the test section of the wind 
tunnel, we placed blocks as roughness for the boundary 
layer. The block size was Sx — 0.06 m, 5y = 0.21m, 
and 8z = 0.11m. The spacing of adjacent blocks was 
Sx — 5y = 0.5 m. 

For the incoming-flow velocities £7* from 2 to 20 ms -1 , 
we measured the streamwise and spanwise velocities at 
x = 12.5 m. The boundary layer was well developed 
there. While the 99% thickness was 0.8 m, the displace- 
ment thickness was 0.2 m f i£ The measurement height was 
z = 0.25-0.35 m in the log-law sublayer. 



III. BURGERS VORTEX 

By using the Burgers vortex, an idealized model for 
vortex tubes, we discuss what information is available 
from a one-dimensional cut of the velocity field. The 
Burgers vortex is an axisymmetric steady circulation in 
a strain field. In cylindrical coordinates, the circulation 
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vortex tube 




uq and strain field (ur,uz) are 



mean stream 



FIG. 1: Sketch of a vortex tube penetrating the (x, y) plane 
at a point (xo, yo)- The inclination is (Go, tpo). The circulation 
velocity is uq. We consider the spanwise velocity v along the 
x axis in the mean stream direction. 
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FIG. 2: Mean profiles in the streamwise («) and spanwise 
(v) velocities for the Burgers vortices with random positions 
(xo,yo) and inclinations (6o,<po)- The u profile is separately 
shown for d x u > and d x u < {u~) at x = 0. The 

position x and velocities are normalized by the radius and 
maximum circulation velocity of the Burgers vortices. The 
dotted line is the v profile of the Burgers vortex for xo = 
yo = 0o = 0, the peak value of which is scaled to that of the 
mean v profile. 
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Here v is the kinematic viscosity and ao (> 0) is a 
constant. The circulation is maximal at R = Rq = 
2.2A{v/ ao) 1 / 2 . We regard Rq as the vortex radius. 

There are other models for vortex tubes, e.g., the 
Lundgren spirals We use the Burgers vortex alone be- 
cause detailed information about individual structures is 
anyway not available from one-dimensional velocity data. 

Suppose that the axis of the Burgers vortex penetrates 
the (x,y) plane at a point (xo,yo) as shown in Fig. [1] 
The x and y axes are in the streamwise and spanwise 
directions. If the inclination of the Burgers vortex is 
(^O)Vo) in spherical coordinates, the streamwise (u) and 
spanwise (v) components of the circulation flow uq along 
the x axis are 



2/oCosflo . . , . 

u{x-x ) = — ue{R), (2a) 

K 

(x-x o )cos0 o , . . , . 

v(x-x ) = ue{R), (2b) 



with 

R 2 = (x - x ) 2 {l - sin 2 9 cos 2 <p ) 

+ yl{\ - sin 2 6» sin 2 v? ) 

-I- 2(x — xo)yo sin 2 #o sin ipo cos <pq. 

Those of the radial inflow ur of the strain field are 



(3) 



(x - x )(l - sin 2 O cos 2 (^o) + yo sin 2 6» sin (p cos /m 

u{x - x ) = u R (R), 

K 

(x - x ) sin 2 6»o sin (^0 cos t^o + J/o(l - sin 2 6» sin 2 ip ) 

v(x - x ) = u R (R). 

K 



(4a) 
(4b) 



When the Burgers vortex is close to the x axis and is not 
heavily inclined, i.e., |yo| ^ Rq and $0 — 0j t ne spanwise 
velocity is dominated by the small-scale circulation flow 
[Eq. (|2b[) ] . The streamwise velocity is dominated by the 
large-scale radial inflow [Eq. (|4a|) ]. When \yo\ 3> Ro or 
80 3> 0, the signal is weak at least for the small-scale 



variation in the spanwise velocity. 

Then, suppose that velocity data are obtained on a 
one-dimensional cut of a flow that consists of vortex tubes 
and the background random fluctuation. The tube posi- 
tion (xo,yo) an( i inclination (6>o,<^o) are supposed to be 
random. This is likely in our experiments, where tur- 
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FIG. 3: Probability density distribution of the absolute 
spanwise-velocity increment \v(x + U/2f s ) — v(x — U/2f 3 )\ 
at Re A = 719, 1304, and 1934. The distribution is verti- 
cally shifted by a factor 10 3 . The increment is normalized by 
{Sv 2 s } 1/2 = ([v(x + U/2f s ) - v(x - U/2f s )] 2 ) 1/2 . The arrows 
indicate the ranges for intense vortical structures, which share 
0.1 and 1% of the total. The dotted line denotes the Gaussian 
distribution. 



FIG. 4: Mean profiles of intense vortical structures for the 
0.1% threshold in the streamwise (u) and spanwise (v) ve- 
locities, (a) Re a = 719. (b) Re A = 1098. The u profile 
is separately shown for d x u > (u + ) and d x u < (u~) 
at x = 0. The position x is normalized by the Kolmogorov 
length. The velocities are normalized by the peak value of the 
v profile. We also show the v profile of the Burgers vortex for 
xo = yo = 6o = by a dotted line. 



bulence was almost isotropic because the measured ra- 
tio (u 2 )/{v 2 ) is not far from unity (Table [J). The vortex 
tubes induce small-scale variations in the spanwise veloc- 
ity. If we consider intense velocity variations above a high 
threshold, their scale and amplitude are close to the ra- 
dius and circulation velocity of intense vortex tubes with 
1 2/o I ^5 Ro an d &o — 0. To demonstrate this, mean profiles 
are calculated for the circulation flows uq of the Burgers 
vortices with random positions (xo,yo) & n d inclinations 
($0) ¥>o)- Their radii Ro and maximum circulation veloc- 
ities Vo = uq(Rq) are set to be the same. We consider 
the Burgers vortices with \d x v\ at x = being above a 
threshold, \d x v\/3 at x — for xq = yo = do = 0. When 
d x v is negative, the sign of the v signal is inverted before 
the averaging. The result is shown in Fig. [5] Despite the 
relatively low threshold, the scale and peak amplitude of 
the mean v profile are still close to those of the v profile 
for xo = yo = So = (dotted line). The extended tails 
are due to the Burgers vortices with |y | ^ Ro or 9o ^S> 0. 



IV. MEAN VELOCITY PROFILE 

Mean profiles of intense vortical structures in the 
streamwise (u) and spanwise (v) velocities are extracted, 



by averaging signals centered at the position where the 
absolute spanwise-velocity increment \v(x + r/2) — v(x — 
r/2)\ is above a threshold i&dlLlld&ll The scale r is the 
sampling interval U/f s . The threshold is such that 0.1% 
or 1% of the increments are used for the averaging (here- 
after, the 0.1% or 1% threshold). These increments com- 
prise the tail of the probability density distribution of all 
the increments as in Fig. O 20 Example of the results are 
shown in Fig. [4}2£ 

The v profile in Fig. Q] is close to the v profile in 
Fig. [2] Hence, the contribution from vortex tubes is 
dominant. The contribution from vortex sheets is not 
dominant. If it were dominant, the v profile should ex- 
hibit some kind of stepJ^ Direct numerical simulations 
at Re A ^5 200 revealed that intense vorticity tends to be 
organized into tubes rather than sheetsi^£&£ ' 21 i 22 This 
tendency appears to exist up to Re> — 2000. Vortex 
sheets might contribute to the extended tails in Fig. 2] 
They are more pronounced than those in Fig. [2j Here 
it should be noted that our discussion is somewhat sim- 
plified because there is no strict division between vortex 
tubes and sheets in real turbulence. 

By fitting the v profile in Fig. Q] around its peaks by 
the v profile of the Burgers vortex for xo = yo = Oo = 
(dotted line), we estimate the radius Ro and maximum 
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TABLE II: Parameters for intense vortical structures: radius Ro, maximum circulation velocity Vo, Reynolds number Reo = 
RoVo/v and small-scale clustering exponent /in- We also list the threshold level tq. 





Units 






Duct flow 










Boundary layer 






1 


2 


3 


4 


5 


6 


7 


s 

O 


n 
■ ) 


10 


11 












For the 0.1% threshold 












TO 


(Sv 2 .) 1 ' 2 
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6.28 
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6.12 


Ro 
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5.21 


5.64 


5.47 


5.87 


6.20 


5.18 


5.24 


5.56 


5.97 


5.66 


6.15 


Vo 


(v 2 ) 1 ' 2 


0.637 


0.595 


0.593 


0.610 


0.629 


0.822 


0.745 
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0.677 
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0.700 
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14.0 
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45.2 
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61.8 


75.1 


86.8 


39.4 


43.9 


53.0 


63.3 


66.2 


78.7 


Mo 




0.867 


0.879 


0.896 


0.919 


0.907 


0.840 


0.855 


0.921 


0.924 


0.995 


0.992 












For the 1% threshold 












TO 


(Sv 2 ) 1 ' 2 


3.66 


3.62 


3.57 


3.55 


3.52 


3.61 


3.64 


3.61 


3.57 


3.57 


3.53 


Ro 


i) 


6.50 


6.61 


6.92 


6.76 


7.03 


6.53 


6.44 


6.92 


7.38 


7.14 


7.19 


Vo 


(v 2 ) 1 ' 2 


0.449 


0.406 


0.395 


0.401 


0.411 


0.607 


0.533 


0.489 


0.466 


0.468 


0.465 


V 


U K 


6.11 


6.83 


7.53 


8.41 


9.17 


5.61 


5.98 


6.65 


7.28 


7.83 


8.53 


Re 




39.7 


45.2 


52.1 


56.9 


64.5 


36.6 


38.5 


46.0 


53.7 


55.9 


61.3 


MO 




1.02 


1.06 


1.07 


1.09 


1.07 


0.884 


0.953 


0.998 


1.03 


1.08 


1.09 



circulation velocity Vq . The measured velocity v m is con- 
sidered to be the true velocity v t averaged over the probe 
size in the streamwise direction, Sx p — 1 mm: 

^ pSx p /2 

v m (x) = - — / v t (x + r)dr. (5) 

0X p J^Sx p /2 

For all the data, the Ro and Vo values are summarized 
in Tabic [TTJ They characterize the scale and intensity 
of vortical structures, even if they are not the Burgers 
vortices. The radius Ro is several times the Kolmogorov 
length 77. The maximum circulation velocity Vo is several 
tenths of the rms velocity fluctuation (v 2 ) 1 / 2 and several 
times the Kolmogorov velocity Similar results were 
obtained from direct numerical simulations^ 4 ^! 7 ^ and 
laboratory experiment o 11 ' 12 ^ 4 ' 15 at the lower Reynolds 
numbers, Re,\ ^ 1300. 

The u profile in Fig. [4] is separated for d x u > 
and d x u < (u~) at x = 0. Since the contamina- 
tion with the w component^ 7 , induces a symmetric posi- 
tive excursion ; 14 i 23 i 24 we decomposed the profiles into 
symmetric and antisymmetric components and show only 
the antisymmetric componentsi 1 ^ The u^ 1 profiles in Fig. 
2] have larger amplitudes than those in Fig. [5] Hence, 
the profiles in Fig. [4] are dominated by the circu- 
lation flows Me of vortex tubes that passed the probe 
with some incidence angles to the mean flow direction, 11 
tan _1 [w/(C7 + u)].@The radial inflow ur of the strain field 
is not discernible, except that the u~ profile has a larger 
amplitude than the u + profilej 14 ' 15 Unlike the Burgers 
vortex, a real vortex tube is not always oriented to the 
stretching direction! 4 ^ 5 .^ 7 -^ 2 ^ 



are dominated by vortex tubes, we expect that the distri- 
bution of intense vortical structures studied here is also 
essentially the distribution of intense vortex tubes. Ex- 
amples of the probability density distribution P{5xo) are 
shown in Figs. [5]and[nP^ 

The probability density distribution has an exponen- 
tial tai l 14 ' 15 that appears linear on the semi-log plot of 
Fig. [5] This exponential law is characteristic of intervals 
for a Poisson process of random and independent events. 
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V. SPATIAL DISTRIBUTION 

The spatial distribution of intense vortical structures 
is studied using the distribution of interval Sxq between 
successive intense velocity increments i 13 i 14 ' 15 i 22 The in- 
tense velocity increment is defined in the same manner 
as for the mean velocity profiles in Sec. HVl Since they 



FIG. 5: Probability density distribution of interval between 
intense vortical structures for the 1% threshold at ReA = 719, 
1304, and 1934. The distribution is normalized by the ampli- 
tude of the exponential tail (dotted line) , and it is vertically 
shifted by a factor 10. The interval is normalized by the 
streamwise correlation length L u . The arrow indicates the 
spanwise correlation length L v . 
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VI. SCALING LAW 
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FIG. 6: Probability density distribution of interval between 
intense vortical structures for the 1% threshold at Ke\ = 719, 
1304, and 1934. The distribution is normalized by the peak 
value, and it is vertically shifted by a factor 10. The dotted 
line indicates the power-law slope from 30?? to 300?;. The 
interval is normalized by the Kolmogorov length r\. The arrow 
indicates the spanwise correlation length L v . 



The large-scale distribution of intense vortical structures 
is random and independent. 

Below the spanwise correlation length L v , the proba- 
bility density is enhanced over that for the exponential 
distribution. 15 Thus, intense vortical structures cluster 
together below the energy-containing scale. In fact, di- 
rect numerical simulations revealed that intense vortex 
tubes lie on borders of energy-containing eddies £ 

Over small intervals, the probability density distribu- 
tion is a power la w 13 ' 22 that appears linear on the log-log 
plot of Fig. H 



P(Sx 



(6) 



Thus, the small-scale clustering of intense vortical struc- 
tures is self-similar and has no characteristic scaled Ta- 
ble |TT] lists the clustering exponent /iq estimated over in- 
tervals from 5xq = 30ry to 3007/. Its value is close to 
unity. 

The exponential law over large intervals and the power 
law over small intervals were also found in laboratory ex- 
periments for regions of low pressure . 26 ' 27 ! 28 ' 29 They are 
associated with vortex tubes, although their radii tend to 
be larger than those of intense vortical structures studied 
here. 29 



Dependence of parameters for intense vortical struc- 
tures on the microscale Reynolds number Re> and on the 
configuration for turbulence production, i.e., duct flow or 
boundary layer, is studied in Fig. [7] Each quantity was 
normalized by its value in the duct flow at Re^ = 1934 
individually for the 0.1% and 1% thresholds. That is, we 
avoid the prefactors that depend on the threshold. When 
the threshold is high, the radius Rq is small, the maxi- 
mum circulation velocity Vq is large, and the clustering 
exponent no is small as in Table HU We focus on scaling 
laws of these quantities. 

The radius Rq scales with the Kolmogorov length rj 
as Rq oc t] [Fig. 0a)]. Thus, intense vortical structures 
remain to be of smallest scales of turbulence. 

The maximum circulation velocity Vq scales with the 
rms velocity fluctuation {v 2 ) 1 / 2 as Vq oc (tj 2 ) 1 / 2 [Fig. 
[TJb)] . Although the rms velocity fluctuation is a charac- 
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FIG. 7: Dependence of parameters for intense vortical struc- 
tures on Re A . (a) Ro/ij. (b) V /{v 2 } 1/2 . (c) Vq/uk- (d) Re . 

1/2 

(e) Reo/Re A . (f) /irj. The open and filled circles respec- 
tively denote the duct flows for the 0.1% and 1% thresholds. 
The upward and downward triangles respectively denote the 
boundary layers for the 0.1% and 1% thresholds. Each quan- 
tity is normalized by its value in the duct flow at Re A = 1934 
individually for the 0.1% and 1% thresholds. 
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teristic of the large-scale flow, vortical structures could 
be formed via shear instability on borders of energy- 
containing eddies,— i 27 ' 28 where a small-scale velocity vari- 
ation could be comparable to the rms velocity fluctua- 
tion. The maximum circulation velocity does not scale 
with the Kolmogorov velocity uk, a characteristic of the 
small-scale flow, as Vq oc uk [Fig. 0(c)]. 

Direct numerical simulations for intense vortex 
tubes^ 7 - at ReA ^ 200 and laboratory experiments for 
intense vortical structure d 11 ! 15 at Re> ^ 1300 derived 
the scalings Rq oc 77 and Vq oc {v 2 ) 1 ^ 2 . We have found 
that these scalings exist up to ReA — 2000, regardless of 
the configuration for turbulence production. 

The scalings of the radius Rq and circulation veloc- 
ity Vo lead to a scaling of the Reynolds number Reo = 
RqVq/v for the intense vortical structures i^ 7 - 

Rc oc Re* /2 if R oc 77 and V oc (v 2 ) 1/2 , (7a) 
Reo = constant if Rq oc rj and Vq oc uk- (7b) 

Our result favors the former scaling [Fig. [71(e)] rather 
than the latter [Fig. [3(d)]. With an increase of ReA, 
intense vortical structures progressively have higher Reo 
and are more unstable^ 7 - Their lifetimes are shorter. It 
is known 3 -^ that the flatness factor ((d x v) 4 ) / ((d x v) 2 } 2 
scales with Re°' 3 . Since ((d x v) ) is dominated by in- 
tense vortical structures, it scales with (v 2 ) 2 /r) 4 . Since 
((d x v) 2 ) 2 is dominated by the background random fluc- 
tuation, it scales with u 4 K jr\ 4 . If the number density of 
intense vortical structures remains the same, we have 
((dxv) 4 ) / ({d x v) 2 ) 2 oc (v 2 ) 2 /u 4 K oc Re 2 ,. The difference 
from the real scaling implies that vortical structures with 
Vo — {v 2 ) 1 / 2 are less numerous at a higher Reynolds num- 
ber Rca, albeit energetically more important. 

The small-scale clustering exponent /io is constant [Fig. 
[3(f)]. A similar result with /io — 1 was obtained from 
laboratory experiments of the Karman flow between two 
rotating disks 2 - at ReA — 400-1600. The small-scale 
clustering of intense vortical structures at high Reynolds 
numbers ReA is independent of the configuration for tur- 
bulence production. 

Lastly, recall that only intense vortical structures are 
considered here. For all vortical structures with vari- 
ous intensities, the scalings Vq oc (v y/ 2 and Re = 

1 /2 

RqVo/v oc Re A are not necessarily expected. For all 
vortex tubes, in fact, direct numerical simulations 3,8 at 
ReA iS 200 derived the scaling Vq oc uk- The devel- 
opment of an experimental method to study all vortical 
structures is desirable. 



VII. CONCLUSION 

The spanwise velocity was measured in duct flows at 
Rca = 719-1934 and in boundary layers at ReA = 332- 
1304 (Table[T|. We used these velocity data to study fea- 
tures of vortical structures, i.e., vortex tubes and sheets. 

We studied the mean velocity profiles of intense vor- 
tical structures (Fig. QJ. The contribution from vortex 
tubes is dominant. Essentially, our results are those for 
vortex tubes. The radius Rq is several times the Kol- 
mogorov length 77. The maximum circulation velocity Vq 
is several tenths of the rms velocity fluctuation (v 2 ) 1 / 2 
and several times the Kolmogorov velocity uk (Table 
ITT]) . There are the scalings Ro oc 77, Vo oc (v 2 ) 1 / 2 , and 
Re = RqVq/v oc Re\ /2 (Fig. EJ). 

We also studied the distribution of interval between in- 
tense vortical structures. Over large intervals, the distri- 
bution obeys an exponential law (Fig. [5]) , which reflects a 
random and independent distribution of intense vortical 
structures. Over small intervals, the distribution obeys a 
power law (Fig. 0, which reflects self-similar clustering 
of intense vortical structures. The clustering exponent is 
constant, ^ - 1 (TableHJand Fig. E}. 

Direct numerical simulations^'^' 7 -' 8 ^'^ and laboratory 
experiment o 11 ' 12 ' 13 ' 14 ! 15 ' 22 derived some of those features. 
We have found that they are independent of the Reynolds 
number and of the configuration for turbulence produc- 
tion, up to ReA — 2000 that exceeds the Reynolds num- 
bers of the prior studies. 

The Reynolds numbers ReA in our study are still lower 
than those of some turbulence, e.g., atmospheric turbu- 
lence at ReA > 10 4 . Such turbulence is expected to 
contain intense vortical structures, because turbulence 
is more intermittent at a higher Reynolds number ReA 
and small-scale intermittency is attributable to intense 
vortical structures. They are expected to have the same 
features as found in our study. These features appear to 
have reached asymptotes at ReA — 2000 (Fig. [7|), regard- 
less of the configuration for turbulence production, and 
hence appear to be universal at high Reynolds numbers 
Re A . 
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